STABILITY AND INSTABILITY OF EXTREME REISSNER-NORDSTROM 
BLACK HOLE SPACETIMES FOR LINEAR SCALAR PERTURBATIONS I 



STEFANOS ARETAKIS* 



Abstract. We study the problem of stability and instability of extreme Reissner-Nordstrom 
spacetimes for linear scalar perturbations. Specifically, we consider solutions to the linear wave 
equation Ogip = on a suitable globally hyperbolic subset of such a spacetime, arising from regular 
initial data prescribed on a Cauchy hypersurface So crossing the future event horizon . We 
obtain boundedness, decay and non-decay results. Our estimates hold up to and including the 
horizon 'H^ . The fundamental new aspect of this problem is the degeneracy of the redshift on 
H'^ ■ Several new analytical features of degenerate horizons are also presented. 

1. Introduction 

Black holes are one of the most celebrated predictions of General Relativity and one of the most 
intriguing objects of Mathematical Physics. A particularly interesting but peculiar example of a 
black hole spacetime is given by the so-called extreme Reissner-Nordstrom metric, which in local 
coordinates {t, r, 6, (p) takes the form 

(1.1) g = ~Ddt^ + ^dr^ +r''gs2, 
where 

gg2 is the standard metric on §^ and M > 0. This spacetime has been the object of considerable 
study in the physics literature (see for instance the recent |38| ) . In this series of papers, we inves- 
tigate the linear stability and instability of extreme Reissner-Nordstrom for scalar perturbations, 
that is to say we shall attempt a more or less complete treatment of the wave equation 

(1.2) agip = o 

on extreme Reissner-Nordstrom exterior backgrounds. This resolves Open Problem 4 (for extreme 
Reissner-Nordstrom) from Section 8 of [5S] . The fundamentally new aspect of this problem is the 
degeneracy of the redshift on the event horizon H"*". Several new analytical features of degenerate 
event horizons are also presented. 

The main results (see Section [4| of the present paper include: 

(1) Local integrated decay of energy, up to and including the event horizon "H^ (Theorem [T|. 

(2) Energy and pointwise uniform boundedness of solutions, up to and including (Theorems 

ii- 

(3) Sharp second order estimates, up to and including "H^ (Theorem^. 

(4) Non-Decay of higher order translation invariant quantities along (Theorem [5|. 

Note that the last result is a statement of instability. In the companion paper [3|, where we shall 
provide the complete picture of the linear stability and instability of extreme Reissner-Nordstrom 
spacetimes, we obtain energy and pointwise decay of solutions up to and including and non- 
decay and blow-up results for higher order derivatives of solutions along H^. Note that the latter 
blow-up estimates are in sharp contrast with the non-extreme case, for which decay holds for all 
higher order derivatives of i(j along . 
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1.1. Preliminaries. Before we discuss in detail our results, let us present the distinguishing prop- 
erties of extreme Reissner-Nordstrom and put this spacetime in the context of previous results. 

1.1.1. Extreme Black Holes. We briefly describe here the geometry of the horizon of extreme 
Reissner-Nordstrom. (For a nice introduction to the relevant notions, we refer the reader to [33]). 



The event horizon corresponds to r = A/, where the (t, r) coordinate system ( 1.1 1 breaks down. 
The coordinate vector field 9t, however, extends to a regular null Killing vector field T on T-L^ . The 
integral curves of T on T-i^ are in fact affinely parametrised: 

(1.3) VtT = 0. 



More generally, if an event horizon admits a Killing tangent vector field T for which (1.31 holds, 



then the horizon is called degenerate and the black hole extreme. In other words, a black hole is 



called extreme if the surface gravity vanishes on 7^+ (see Section 2.3 ). Under suitable circumstances, 
the notion of extreme black holes can in fact be defined even in case the spacetime does not admit 
a Killing vector field (see [20] )• There has also been rapid progress as regards the problem of 
uniqueness of extreme black hole spacetimes. We refer the reader to [16] for the classification of 
static electro-vacuum solutions with degenerate components. 

The extreme Reissner-Nordstrom corresponds to the M — e subfamily of the two parameter 
Reissner-Nordstrom family with parameters mass M > and charge e > 0. It sits between the 
non-extreme black hole case e < M and the so-called naked singularity case M < e. Note that the 
physical relevance of the black hole notion rests in the expectation that black holes are "stable" 
objects in the context of the dynamics of the Cauchy problem for the Einstein equations. On the 
other hand, the so-called weak cosmic censorship conjecture suggests that naked singularities are 
dynamically unstable (see the discussion in |26)V That is to say, extreme black holes are expected to 
have both stable and unstable properties; this makes their analysis very interesting and challenging. 

1.1.2. Linear Scalar Perturbations. The first step in understanding the non- linear stability (or in- 



stability) of a spacetime is by considering the wave equation (1.2) (we refer the reader to [TT] for 
the details of the proof of the stability of Minkowski space). This is precisely the motivation of 
the present paper. Indeed, to show stability one would have to prove that solutions of the wave 
equation decay sufficiently fast. For potential future applications all methods should ideally be 
robust and the resulting estimates quantitative. Robust means that the methods still apply when 
the background metric is replaced by a nearby metric, and quantitative means that any estimate 
for must be in terms of uniform constants and (weighted Sobolev) norms of the initial data. Note 
also that it is essential to obtain non-degenerate estimates for tp on H"*" and to consider initial data 
that do not vanish on the horizon. As we shall see, the issues at the horizon turn out to be the 
most challenging part in understanding the evolution of waves on extreme Reissner-Nordstrom. 



1.1.3. Previous Results for Waves on Non-Extreme Black Holes. The wave equation (1.2) on black 
hole spacetimes has been studied for a long time beginning with the pioneering work of Regge and 
Wheeler [JS] for Schwarzschild. Subsequently, a series of heuristic and numerical arguments were 
put forth for obtaining decay results for ip (see [6l|43]). However, the first complete quantitative 
result (uniform boundedness) was obtained only in 1989 by Kay and Wald [56], extending the 
restricted result of |51j . Note that the proof of 56 heavily depends on the exact symmetries of the 
Schwarzschild spacetime. 

During the last decade, the wave equation on black hole spacetimes has become a very active 
area in mathematical physics. As regards the Schwarzschild spacetime, "X estimates" providing 



local integrated energy decay (see Section 1.2.1 below) were derived in [7l|8l|2T]. Note that [21] 



introduced a vector field estimate which captures in a stable manner the so-called redshift effect, 
which allowed the authors to obtain quantitative pointwise estimates on the horizon H.^ . Refine- 
ments for Schwarzschild were achieved in and fS^. Similar estimates to (Tj were derived in [S] 
for the whole parameter range of Reissner-Nordstrom including the extreme case. However, these 
estimates degenerate on H"*" and require the initial data to be supported away from H.^ . 

The first boundedness result for solutions of the wave equation on slowly rotating Kerr (|a| <C M) 
spacetimes was proved in [23j and decay results were derived in ^ ,25t i50j . Decay results for general 
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subextreme Kerr spacetimes (|a| < M) are proven in 129 . Two new methods were presented recently 
for obtaining sharp decay of energy flux and pointwise decay on black hole spacetimes; see [?7ll51) . 
For results on the coupled wave equation see [5^. For other results see [inilUlISS]. For an exhaustive 
list of references, see [25] , 

Note that all previous arguments for obtaining boundedness and decay results on non-extreme 



black hole spacetimes near the horizon would break down in our case (see Sections 2.3 10.1 1. The 
reason for this is precisely the degeneracy of the redshift on H+. 

1.2. Overview of Results and Techniques. We use the robust vector field method (see Section 
[5]). Our methods at various points rely on the spherical symmetry but not on other unstable 
properties of the spacetimes (such as the complete integrability of the geodesic flow, the separability 
of the wave equation, etc.) 

1.2.1. Zeroth Order Morawetz and X Estimates. Our analysis begins with local spacetime esti- 
mates. We refer to local spacetime estimates controlling the derivatives of i/j as "X estimates" and 
■0 itself as "zeroth order Morawetz estimate" . Both these types of estimates have a long history 
(see [25]) beginning with the seminal work of Morawetz |41j for the wave equation on Minkowski 
spacetime. They arise from the spacetime term of energy currents associated to a vector field X 
(see Section [5] for the definition of energy currents). For Schwarzschild, such estimates appeared in 
OIHIIIO, 21, 23 and for Rcissner-Nordstrom in 9]. The biggest difficulty in deriving an X estimate 
for black hole spacetimes has to do with the trapping effect. Indeed, from a continuity argument 
one can infer the existence of null geodesies which neither cross nor terminate at In our 
case, a class of such geodesies lie in a hypersurface of constant radius (see Section |2.2[ ) known as 
the photon sphere. From the analytical point of view, trapping affects the derivatives tangential 
to the photon sphere and any non-degenerate spacetime estimate must lose (tangential) derivatives 
(i.e. must require high regularity for 

In this paper, we first (making minimal use of the spherical decomposition) derive a zeroth order 
Morawetz estimate for '4> which docs not degenerate at the photon sphere. For the case I > 1 (where 
I is related to the eigenvalues of the spherical Laplacian, see Section [7| we use a suitable current 
which captures the trapping in the extreme case. As regards the zeroth spherical harmonics {I ~ 0), 
we present a method which is robust and uses only geometric properties of the domain of outer 
communications. Our argument (for I = 0) applies for a wider class of black hole spacetimes and, 
in particular, it applies for Schwarzschild. Note that no unphysical conditions are imposed on the 
initial data which, in particular, are not required to be compactly supported or supported away 
from H^. Once this Morawetz estimate is established, we then show how to derive a degenerate 
(at the photon sphere) X estimate which does not require higher regularity and a non-degenerate 
X estimate (for which we need, however, to commute with the Killing vector field T). These 



estimates, however, degenerate on this degeneracy will be eliminated later (see Section 1.2.4 1. 
See Theorem [T] of Section |4l 

1.2.2. Uniform Boundedness of Non-Degenerate Energy. The vector field T ^ dt is causal and 
Killing and the energy fiux of the current is non-negative definite (and bounded) but degenerates 
on the horizon (see Section [8|. Moreover, in view of the lack of redshift along H+, the divergence 
of the energy current associated to the redshift vector field N, first introduced in [21], is not 
positive definite near (see Section 10). For this reason we appropriately modify so the new 
hulk term is non-negative definite near H"*". Note that the bulk term is not positive far away from 
T-i^ and so to control these terms we use the X and zeroth order Morawetz estimates. The arising 
boundary terms can be bounded using Hardy-like inequalities. It is important here to mention that 
a Hardy inequality (in the first form presented in Section^^ allows us to hound the local norm 
of ij; on hypersurfaces crossing T-L^ using the (conserved) degenerate energy ofT. See Theorem [2] 
of Section m 

1.2.3. Non-Decay along . Wc next show that the degeneracy of redshift gives rise to a conser- 
vation law along for the zeroth spherical harmonics. This result implies that for generic waves 
a higher order translation invariant geometric quantity does not decay along 'H'^ (see Theorem [s] of 
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Section [4]). This low frequency obstruction will be crucial for obtaining the definitive statement of 
instability. 

1.2.4. Local Integrated Energy Decay. We next return to the problem of retrieving the derivative 
transversal to "H"*" in the X estimate in a neighbourhood of H'^; see Theoremjsjof Sectionji] We first 
show that on top of the above low frequency obstruction comes another new feature of degenerate 
event horizons. Indeed, to obtain non-degenerate spacetime estimates, one needs to require higher 
regularity for tp and commute with the vector field transversal to T-L^ . This shows that T-L^ exhibits 



phenomena characteristic of trapping (see also the discussion in Section 1.3.2 1. Then by using 
appropriate modifications of the redshift current and Hardy inequalities along T-i^ we obtain the 
sharpest possible result. See Section [l2j Note that although (an appropriate modification of) the 
redshift current can be used as a multiplier for all angular frequencies, the redshift vector field N 
can only be used as a commutatoi]^ for ip supported on the frequencies ^ > 1. These results will be 
further investigated in [3]. 

1.2.5. Pointwise Uniform Boundedness. Using the above higher order energy estimates and appro- 
priate Sobolev inequalities we finally obtain uniform pointwise boundedness of solutions up to and 
including T-i^ (see Theorem |4] of Section |4]). We note that in [3] we show that the argument of Kay 
and Wald |56| cannot be applied in the extreme case for obtaining similar boundedness results, i.e. 
for generic if), there does not exist a Cauchy hypersurface S crossing 7^+ and a solution ij; such that 

Til) — Tp 

in the causal future of S. 

1.3. Remarks on the Analysis of Extreme Black Holes. We conclude this introductory sec- 
tion by discussing several new features of degenerate event horizons. 

1.3.1. Dispersion vs Redshift. In |25) . it was shown that for a wide variety of non-extreme black 
holes, the redshift on suffices to yield uniform boundedness (up to and including the horizon) of 
waves ip without any need of understanding the dispersion properties of ^p. However, in the extreme 
case, the degeneracy of the redshift makes the understanding of the dispersion of ip essential even 
for the problem of boundedness. In particular, one has to derive spacetime integral estimates for 
-0 and its derivatives. Moreover, we show that dispersion can be completely decoupled from the 
redshift effect. 



1.3.2. Trapping Effect on H'^ . According to the results of Sections 10.2 and 12 in order to obtain 
estimates in neighbourhoods of TL^ one must require higher regularity for ^p and commute with 
the vector field transversal to (which is not Killing). This loss of a derivative is characteristic of 
trapping. Geometrically, this is related to the fact that the null generators of H'^ viewed as integral 
curves of the Killing vector field T are affinely parametrised. 

The trapping properties of the photon sphere have a different analytical flavour. Indeed, in 
order to obtain estimates in regions which include the photon sphere, one needs to commute 
with either T or the generators of the Lie algebra so (3) (note that all these vector fields are Killing). 
Only the high angular frequencies are trapped on the photon sphere (and for the low frequencies 
no commutation is required) while all the angular frequencies are trapped (in the above sense) on 

n+. 

2. Geometry of Extreme Reissner-Nordstrom Spacetime 

The unique family of spherically symmetric asymptotically flat solutions of the coupled Einstein- 
Maxwell equations is the two parameter Reissner-Nordstrom family of 4-dimensional Lorentzian 
manifolds (AfM,e, gM.e) where the parameters M and e are called mass and (electromagnetic) charge, 
respectively. 

The Reissner-Nordstrom metric was first written in local coordinates (t, r, 6, (p) in 1916 46 and 
1918 [12]. In these coordinates, g = gM,e = -Ddt^ + j^dr^ + r^g^2 , where D = D (r) = 1 - ^ + 



^The redshift vector field was used as a commutator for the first time in |23 
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and gg2 is the standard metric on S"^. The extreme case corresponds to M = |e|. From now on we 
restrict our attention to the extreme case. 

Clearly, S0(3) acts by isometry on these spacetimes. We will refer to the SO(3)-orbits as (symme- 
try) spheres. The coordinate r is defined intrinsically such that the area of the spheres of symmetry 
is 47rr^ (and thus should be thought of as a purely geometric function of the spacetime). 

In view of the coordinate singularity a,t r — M (note that M is the double root of D in the 
extreme case), we introduce the so-called tortoise coordinate r* given by ^^g^^'' — j^- Note that in 
the extreme case r* is inverse linear (instead of logarithmic in the non-extreme case). The metric 
with respect to the system {t, r*) then becomes g = —Ddt^ + D (dr*) + r'^g^^. 

The way to extend the metric beyond r = M is by considering the ingoing Eddington-Finkelstein 
coordinates (v, r) where v = t + r* . In these coordinates the metric is given by 

(2.1) g = ~Ddv^ + 2dvdr + r^g^2 , 

where D = D (r) ~ (l — and g^2 is the standard metric on S^. The radial curves v — c, 
where c is a constant, are the ingoing radial null geodesies. This means that the null coordinate 
vector field 9,. differentiates with respect to r on these null hypersurfaces. This geometric property 
of dr makes this vector field very useful for understanding the behaviour of solutions to the wave 
equation close to "H^. 

The Penrose diagram of the spacetime Af covered by this coordinate system for w e M, r G M+ is 




We will refer to the hypersurface r = M as the event horizon (and denote it by H~^) and the 
region r < M as the black hole region. The region where M < r corresponds to the domain of 
outer communications. 

In view of the existence of the timelike curvature singularity {r = 0} 'inside' the black hole 
(thought of here as a singular boundary of the black hole region) and its unstable behaviour, we are 
only interested in studying the wave equation in the domain of outer communications including 
the horizon H"*". Note that the study of the horizon is of fundamental importance since any 
attempt to prove the nonlinear stability or instability of the exterior of black holes must come to 
terms with the structure of the horizon. 

We consider a connected asymptotically fiat S0(3)-invariant spacelike hypersurface So in Af 
terminating at i° with boundary such that dT,Q — T^qC] . We also assume that if n is its future 
directed unit normal and T = dy then there exist positive constants Ci < C2 such that 

Ci < -g (n, n) < C2, Ci < -g {n, T) < C2. 
Let M be the domain of dependence of Sq- Then, using the coordinate system (w,r) we have 
(2.2) M = ((-00, +00) X [M , -1-00 ) X §2) n J+ (So) , 

where J+ (Eq) is the causal future of Sq (which by our convention includes Eo)- Note that is a 
manifold with stratified (piecewise smooth) boundary dM = {T-L^ n 7W) U Eo- 
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~cf.'_ 
I 



Another coordinate system that partially covers M. is the null system (u, v) where u — t — r* ,v ^ 
t + r* and with respect to which the metric is g — —Ddudv + r^g^^. The hypersurfaces v = c and 
u — c are null and thus this system is useful for applying the method of characteristics. 

2.1. The Foliations and S,-. We consider the foliation St- = (pr{^o), where ipr is the flow of 
T = dy. Of course, since T is Killing, the hypersurfaces S,- are all isometric to Sq. 



i+ 




We define the region 

71(0, t) = Uo<f<TSf . 

On Tir we have an induced Lie propagated coordinate system such that p G [M,+oo) and 

cj e These coordinates are defined such that if Q G S,- and Q = {vQ,rQ,ujQ) then p = rg and 
id — ujQ. Our assumption for the normal (and thus for ris^) implies that there exists a bounded 
function gi such that 

(2.3) dp^gidi,+dr. 

This defines a coordinate system since [dp, de] = [dp, 9^] = 0. Moreover, the volume form of T,r is 

(2.4) dgs, = Vp^dpduj, 
where ^ is a positive bounded function. 




In the companion paper [3j we shall make use of another foliation denoted by E^. whose leaves 
terminate at 1+ and thus "follow" the waves to the future. 
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i+ 




One can similarly define an induced coordinate system on E,-. Note that only local elliptic 
estimates are to be applied on Et-. 

2.2. The Photon Sphere and Trapping Effect. One can easily see that there exist orbiting 
future directed null geodesies, i.e. null geodesies that neither cross the horizon nor meet null 
infinity 1+ . A class of such geodesies 7 are of the form 

7 :M ^ X 

^->7(r)- (t(T),Q,|,0(r)). 

The conditions V.^7 = and g (7, 7) = imply that Q — 2M, which is the radius of the so called 
photon sphere. The t, cf) depend linearly on r. 



In fact, from any point there is a codimension-one subset of future directed null directions whose 
corresponding geodesies approach the photon sphere to the future. The existence of this "sphere" 
(which is in fact a 3-dimensional timelike hypersurface) implies that the energy of some photons is 
not scattered to null infinity or the black hole region. This is the so called trapping effect. As we 
shall see, this effect forces us to require higher regularity for the waves in order to achieve decay 
results. 

2.3. The Redshift Effect and Surface Gravity of The Killing vector field dy becomes 
null on the event horizon "H^ and is also tangent to it and thus T-L^ is a Killing horizon. In general, 
if there exists a Killing vector field V which is normal to a null hupersurface then 

(2.5) VvV^kV 

on the hypersurface. Since V is Killing, the function k is constant along the integral curves of V . 
This can be seen by taking the pushforward of ( |2.5[ ) via the fiow of V and noting that since the flow 
of V consists of isometries, the pushforward of the Levi-Civita connection is the same connection. 
The quantity k is called the surface gravit"^ of the null hypersurfac^ In the Reissner-Nordstrom 
family, the surface gravities of the two horizons {r = r_} and {r = r+j are given by 

(2.6) ^^_r^~r^ _ldD{r) 



2r| 2 dr 



where D is given in Section [2] (and r+,r_ are the roots of D). Note that in extreme Reissner- 
Nordstrom spacetime we have D{r) = (l — ^) and so r^ — = M which implies that the 
surface gravity vanishes. A horizon whose surface gravity vanishes is called degenerate. 



■^This plays a significant role in black hole "thermodynamics" (see also 1551 and 1441 ). 

''Note that in Riemannian geometry, any Killing vector field that satisfies ||2.5^ must have k = 0. 
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Physically, the surface gravity is related to the so-called redshift effect that is observed along (and 
close to) H^. According to this effect, the wavelength of radiation close to Ti^ becomes longer as 
V increases and thus the radiation gets less energetic. This effect has a long history in the heuristic 
analysis of waves but only in the last decade has it been used mathematically. For example, Price's 
law (see |20j ) and the stability and instability of Cauchy horizons in an appropriate setting (see 
[IH]) were proved using, in particular, this effect. (Note that for the latter, one also needs to use 
the dual blueshift effect which is present in the interior of black holes.) 



3. The Cauchy Problem for the Wave Equation 



We consider solutions of the Cauchy problem of the wave equation (1.2) with initial data 

(3.1) V-lso = V-o e Ht, (So) , nsoVlso = V'l ^ ^^it' (^o) , 

where the hypersurface Eg is as defined in Section [2] and ti^q denotes the future unit normal of 
Sq. In view of the global hyperbolicity of M, there exists a unique solution to the above equation. 
Moreover, as long as fc > 1, we have that for any spacelike hypersurface S 

In this paper we will be interested in the case where k > 2. Moreover, we assume that 

(3.2) lim r^^{x) = 0. 

For simplicity, from now on, when we say "for all solutions ip of the wave equation" we will 
assume that tp satisfies the above conditions. Note that for obtaining sharp decay results we 
will have to consider even higher regularity for ijj. 



4. The Main Theorems 

We consider the Cauchy problem for the wave equation (see Section [s]) on the extreme Reissner- 
Nordstrom spacetime. This spacetime is partially covered by the coordinate systems {t,r), {t,r*), 
{v,r) and {u,v) described in Section [2] Recall that Af is a positive parameter and D — D{r) = 
(l — ^) . Recall also that the horizon Ti^ is located at {r = M} and the photon sphere at 
{r = 2M}. 

We denote T ~ dy ^ dt, where dy corresponds to the system (f,r) and dt corresponds to (t, r). 
From now on, d^, dr are the coordinate vector fields corresponding to (w, r), unless otherwise 
stated. Note that dr- = 9t, = T on Ti.'^ and therefore it is not transversal to whereas dr is 
transversal to Ti,^ . 

The foliation Yir is defined in Section 2.1 and the current associated to the vector field V is 



defined in Section |5.1| Note that every time we use such a current we refer to I^ as a multiplier 
vectorfield. For reference, we mention that 

jjiV'in^ - im' + (^1 - ^Jidr^^r + , 

which degenerates on whereas 



2 



which does not degenerate on H^. The Fourier decomposition of t/j on S^(r) is discussed in Section 
[7] where it is also defined what it means for a function to be supported on a given range of angular 
frequencies. Note that all the integrals are considered with respect to the volume form. The initial 
data are assumed to be as in Section |3] and sufficiently regular such that the right hand side of the 
estimates below are all finite. Then we have the following 

Theorem 1. (Morawetz and X Estimates) Let 6 > 0. There exists a constant Cs > which 
depends on M , 5 and Sg such that for all solutions tp of the wave equation the following estimates 
hold 
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(1) Non-Degenerate Zeroth Order Morawetz Estimate: 



L 



K(0,r) ^ JT,o 

(2) X Estimate with Degeneracy at and Photon Sphere: 

'(r-2A-/)2./D 



/ 



where x^ml 5m]('') *s the indicator function of the interval [ 
(3) X Estimate with Degeneracy at : 



3M 5M 1 
2 ' 2 J 



k(o,t)\ ^ ^ r I I- 



For proving statement (1), we first derive in Section 9.7 an estimate which degenerates on "H^ 



Note that the estimate of Section 9.7 requires applying only the Killing vectorfield T as a multiplier. 
Later in Section [M] we construct a timelike multiplier and we use an appropriate version of Hardy's 
inequality to eliminate this degeneracy (see Theorem [s]). 

Note that estimate (2) degenerates on the photon sphere with respect to all derivatives except 
precisely dr- ■ Note moreover, that the same estimate degenerates on "H^ with respect to all deriva- 
tives. As we shall see, unlike the subextreme case, the degeneracy on with respect to the 
transversal derivative dr can not be eliminated even if we apply a timelike multiplier (see also 
Theorem [2]). 

In Section [9.6| we commute with the vectorfield T to retrieve the tangential to the photon sphere 
derivatives thus obtaining statement (3). 

Theorem 2. (Uniform Boundedness of N on- Degenerate Energy) There exists ro such that 
M < rQ < 2M and a constant C > which depends on M, and Eq such that if A ^ 7?.(0,t) n 
{M < r < To} then for all solutions ip of the wave equation we have 



-' ^ -'So 



The above theorem is proved in Section [TT] using a novel redshift current constructed in Section 
10 and a new version of Hardy's inequality (see Section [6]). Note that this theorem shows that 
we can not eliminate the degeneracy of dr in spacetime neighbourhoods of T-L^ even if we apply 
timelike multipliers. To eliminate this degeneracy we need to reveal a new feature of degenerate 
horizons captured in the next theorem. 

Theorem 3. (Trapping Effect on the Event Horizon H^) Let A be the spacetime region 
defined in Theorem^ and 5 > 0. Then there exists a constant C > (and Cg > 0) which depends 
on M and Eg (and S) such that for all solutions ip with vanishing spherical mean (i.e. for all tp 
supported on the angular frequencies I > 1), the following hold 

(1) Sharp Second Order Estimates: 

{d^drlPf + {drdrlPf + 1^9,^1' + / {d^dri'f + Xl 1^5,^1' 



< C / [i;]n^^^ + C / J^" [TV-ln^^ + C / J^" , 

"'So "'So "'SoH.A 

where xi = i/';/' is supported on I = 1 and xi = 1 i/V' is supported on I >2. 
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(2) Local Integrated Energy Decay: 

Statement (1) of Theorem |3] is proved in Section 12 where we construct a new current and 
derive appropriate Hardy inequahties along T-L^ . Note that the restriction on frequencies / > 1 is 
required in view of a new low- frequency phenomenon of degenerate horizons (see also Theorem [5|. 
In particular, we note that in [S, we show that there is no constant C such that statement (1) holds 
for all solutions ij: of the wave equation. Therefore, the assumption on the frequency range is sharp. 

In statement (2) we have eliminated the degeneracy of dr in neighbourhoods of at the 
expense, however, of commuting the wave equation with dr and thus requiring higher regularity for 
ip. This allows us to conclude that the event horizon T-i^ exhibits trapping. 

Theorem 4. (Pointwise Boundedness) There exists a constant C which depends on AI and Sq 
such that for all solutions ip of the wave equation we have 



everywhere in TZ, where 



1^1 < cJe^, 



./En -/So 



This theorem is proved in Section [M] using Theorem [2] and Sobolev inequalities. 

Theorem 5. (Non-Decay) For generic initial data which give rise to solutions ijj of the wave 
equation the quantity 

does not decay along Ti,^ . 

This theorem arises from a new low-frequency phenomenon of degenerate horizons. In particular, 
in Section [12.1| we show that an appropriate expression which corresponds to the spherical mean of 
■0 is conserved along T-i^ under the evolution. This low-frequency "instability" will be extensively 
investigated in the companion papers O 0] . 

5. The Vector Field Method 

For understanding the evolution of waves we will use the so-called vector field method. This is 
a geometric and robust method and involves mainly estimates. The main idea is to construct 
appropriate (0,1) currents and use Stokes' theorem (see Appendix [b|) in appropriate regions. For a 
nice recent exposition see [1 . We briefly recall here the method to set notation. 

Given a (0,1) current we have the continuity equation 

(5.1) / P,<^ = / P,<^ + [ PXn+ + f V^P, 

JSo JT.^ J-H+{0,t) JtZ{0,t) 

where all the integrals are with respect to the induced volume form and the unit normals n^^ are 
future directed. All the integrals are to be considered with respect to the induced volume form and 
thus we omit writing the measure. Our choice is — dy — T (and thus the volume element on 



y.^ is chosen so (5.1) holds) 



5.1. The Compatible Currents J, K and the Current £. We usually consider currents P^ 
that depend on the geometry of {M,g) and are such that both P^ and V^P^ depend only on the 
1-jet of ip. This can be achieved by using the wave equation to make all second order derivatives 
disappear. There is a general method for producing such currents using the energy momentum 
tensor T. Indeed, the Lagrangian structure of the wave equation gives us the following energy 
momentum tensor 

(5.2) T^, [^] = S^VS.V - Ig^^.d^'i^dai^, 
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which is a symmetric divergence free (0,2) tensor. We wiU in fact consider this tensor for general 
functions tp : A4 ^ R in which case we have the identity 

(5.3) DivT = (Og^) dip. 

Since T is a (0,2) tensor we need to contract it with vector fields of A^. It is here where the 
geometry of Ai makes its appearance. Given a vector field V we define the current by 

(5.4) J^[^]^T,MV''. 

We say that we use the vector field ^ as a 77iM/izp/ie7|^ if we apply ( |5.1| ) for the current . The 

divergence of this current is Div(J) = Div (T) y + T (V^) , where (VVY^ = {g'^'VkVy = {V'VY . 
If ip is a, solution of the wave equation then Div (T) — and, therefore, V'^ is an expression of 
the 1-jet of tp. Given a vector field V the scalar current is defined by 

(5.5) if ^ m = T [V.] (vv) = T,, [ip] {v^vy . 

Note that from the symmetry of T we have [ip] — Tf^^, [ip] TTy" , where iTy^ ~ [Cvg]^'^ is the 
deformation tensor of V . Clearly if ip satisfies the wave equation then if ^["0] = J^['0]. Thus if we 
use Killing vector fields as multipliers then the divergence vanishes and so we obtain a conservation 
law. This is partly the content of a deep theorem of Noethei]^ In general we define the scalar 
current by 

(5.6) [0] = Div (T) = {UgiP) diPiV)=^ (Dg V) (V^) ■ 

5.2. The Hyperbolicity of the Wave Equation. The hyperbolicity of the wave equation is 
captured by the following proposition 

Proposition 5.2.1. Let Vi, V2 be two future directed timelike vectors. Then the quadratic expression 
r(Vi, V2) is positive definite in dtp. By continuity, if one of these vectors is null then T{Vi^ V2) is 
non-negative definite in dip. 

Proof. Gonsider a point p ^ M and the normal coordinates around this point and that without 
loss of generality V2 — (1, 0, 0, 0). Then the proposition is an application of the Cauchy-Schwarz 
inequality. □ 

An important application of the vector field method and the hyperbolicity of the wave equation 
is the domain of dependence property. 

As we shall see, the exact dependence of T on the derivatives of ip will be crucial later. For a 
general computation see Appendix |A.2[ 



6. Hardy Inequalities 

In this section we establish three Hardy inequalities which as we shall see will be very crucial 
for obtaining sharp estimates. They do not require ip to satisfy the wave equation and so we only 
assume that ip satifies the regularity assumptions described in Section [s] and (3.2). From now on, 



St- is the foliation introduced in Section 2.1 (however, these inequalities hold also for the foliation 

tr). 

Proposition 6.0.2. (First Hardy Inequality) For all functions ip which satisfy the regularity 
assumptions of Section^we have 

1 



r2 



where the constant C depends only on M and Sq. 



the name comes from the fact that the tensor T is multiplied by V. 

'According to this theorem, any continuous family of isometrics gives rise to a conservation law. 
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Proof. Consider the induced coordinate system {p,uj) introduced in Section 2.1 We use the 1- 
dimensional identity 

/•+00 r+oo 

(6.1) / {dph)^^dp^[h4,XZi-'^ h^{dp^)dp 

Jm Jm 
with h — r — M. In view of the assumptions on ip, the boundary terms vanish. Thus, Cauchy- 
Schwarz gives 

i;^dp = -2 [ {r^M)^{dp^)dp<2[ ( ^?dp\ If {r - Mf{dpi})^dp] . 

{r>M} J{r>M} \-J {r>M} J \J {r>M} J 

Therefore, 

i'^dp<Aj {r- Mfidpipfdp. 

{r>M} J{r>M} 



Integrating this inequality over §^ gives us 

1 



il?p^dpduj <'i i / D{dpi!fp^dpduj. 

§2 J{r>M} P Js^ J{r>M} 



The result follows from the fact that each S,- is diffeomorphic to §^ x [M, +oo), from (2.3) and the 



boundedness of the factor V in the volume form (2.4) of E^- D 



The importance of the above inequality lies in the weights. The weight of (dptp)'^ vanishes to 
second order on but does not degenerate at infinitjl^ whereas the weight of ip degenerate at 
infinity but not at r = M. Similarly, one might derive estimates for the non-extreme cas^ We also 
mention that the right hand side is bounded by the (conserved) flux of T through E^. (see Section 

Proposition 6.0.3. (Second Hardy Inequality) Let tq E {M,2AI). Then for all functions 
which satisfy the regularity assumptions of Section and any positive number e we have 

^^<e f (9„^)2 + (9.V)' +Ce f 

where the constant depends on M, e, and Eg. 



Proof. We use as before the 1-dimensional identity (6.1 ) with h = r — r^. Then 



(ro - M)V''(M) = / ° ^2 ^ 2g^{dp^)dp <e [" [dp^jfdp + [" + i'^dp, 
Jm jm Jm V ^ 



for any e > 0. By integrating over § , using (2.3) and noting that the p factor that appears in the 



volume form of E,- H {r < rg} is bounded we obtain the required result. □ 

The previous two inequalities concern the hypersurfaces E,- that cross . The next estimate 
concerns spacetime neighbourhoods of "H^. 

Proposition 6.0.4. (Third Hardy Inequality) Let ro,ri be such that M < tq < ri. We define 
the regions A = TZ{0, r) n {M < r < rg} , B = 7?.(0, r) n {rg < r < ri} . Then for all functions ip 
which satisfy the regularity assumptions of Section we have 

ij^<c [ +c [ D[{d,ijf + {dr4>)\ 

A Jb Jaub 

where the constant C depends on M , rp, ri and Eq. 



''Note that if we replace r — M with another function h then we will not be able to make this weight degenerate 
fast enough without obtaining non-trinial boundary terms. 

^These estimates turn out to be stronger since the weight of ^ may diverge at r = M in an integrable manner. 
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Proof. We use again (6.1) with h such that h ~ 2{p — M) in [A/, tq] and h{ri) = 0. Then 

pro pro pri nri 

2/ ^^dp< ^^dp+ (i-(dpg))^^dp+ g^dp^fdp. 

Jm Jm Jro JM 

Thus, integrating over we obtain 

^^<C f D[{dy^f + {dr^f] +c f 

for aU T > 0, where the constant C depends on M, tq, ri and Eq- Therefore, integrating over 
T G [0, r] and using the coarea formula 



'0 KJSf 

completes the proof of the proposition. 



f]dT 



f 



□ 



7. Elliptic Theory on §2(r),r > 

In view of the symmetries of the spacetime, it is important to understand the behaviour of 
functions on the orbits of the action of SO (3), which are isometric to §^(r) for r > 0. Recall that 
the eigenvalues of the spherical Laplacian /f\ are equal to ^^^5^,^ G N. The dimension of the 
eigenspaces is equal to 2^ + 1 and the corresponding eigenvectors are denoted by y™'',— Z < 
m < I and called spherical harmonics. We have L^(§^(r)) = ©oo-E' and, therefore, any function 
-0 G L^(S^(r)) can be written as 



00 / 



(7.1) 



i=0 m=-l 



The right hand side converges to ip in of the spheres and under stronger regularity assumptions 
the convergence is pointwise. Let us denote by ipi the projection of %l) onto E\ i.e. 



Since each eigenspace is finite dimensional (and so complete) and so closed we have L^(§^(r)) = 

E'- ® [E^) . Therefore, ip can be written uniquely as 1/' = V-'o + ''/'>ii where G i-^^") = ®i>iE^ ■ 
We also have the following 

Proposition 7.0.5. (Poincare inequality) If ip € (§^(r)) and tpi = for all I < L — 1 for 

some finite natural number L then we have 



L[L + 1) 



< / lyvr 



and equality holds if and only if ipi — for all I ^ L. 
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Proof. We have ip — J2i>l Therefore, 

^ L(£ + l) 




where we have used the orthogonahty of distinct eigenspaces. □ 
Returning to our 4-dimensional problem, if tp is sufficiently regular in TZ then its restriction at 



each sphere can be written as in ( |7.1[ ), where V'm,; = tpm,i{'f, r) and F™'' = Y™-'\6, </>). We will not 
worry about the convergence of the series since we may assume that ip is sufficiently regulai]^ First 



observe that for each summand in (|7.1|) we have 

in I 1^ \ 



where 5* is an operator on the quotient A^/S0(3). Therefore, if ■(/; satisfies the wave equation 
then in view of the linear independence of y™>''s the terms i/'m,; satisfy S'V'm.i = ^^Tr^V'm,; and, 
therefore, each summand also satisfies the wave equation. This implies that if ?A„i^; — initially then 
■^m,; = everywhere. From now on, we will say that the wave i/; is supported on the angular 
frequencies / > L if t/j^ = 0, i = 0, 1 initially (and thus everywhere). Similarly, we 

will also say that t/; is supported on the angular frequency I ^ L ii ip E . 

Another important observation is that the operators dy and dr are endomorphisms of for all 
I. Indeed, dy commutes with /fs and \i E E^ then 

and thus € ■ Therefore, if 4' is supported on frequencies I > L then the same holds for drip. 

8. The Vector Field T 

One can easily see that = dt = df in the intersection of the corresponding coordinate systems. 
Here dt corresponds to the coordinate basis vector field of either (<, r) or (i, r*). Recall that the 
region M. where we want to understand the behaviour of waves is covered by the system (v,r,9, (jj). 



Therefore we define T = dy. It can be easily seen from (2.11 that T is Killing vectorfield and 
timelike everywhere^ except on the horizon where it is null. 

8.1. Uniform Boundedness of Degenerate Energy. Recall that — T^i^tt^'', where tt^"^ is 
the deformation tensor of T. Since T is Killing vectorfield, its deformation tensor is zero and so 
— 0. Therefore, the divergence identity in the region Tl{0, r) gives us the following conservation 

law 

(8.1) / j^mn'^^+ f j^mn^^^^ f JJ[^K„. 

J'H+ Jso 

Since T is null on the horizon H'^ we have that J^nt^+ > 0. More presicely, since n!^+ = T, from 
(A.l) (See Appendix A. 2) we have Jj^nt^^ — (dyip)'^, thus proving the following proposition: 



Proposition 8.1.1. For all solutions ip of the wave equation we have 
(8.2) / jjM^i < / J^i^Wi,. 



Indeed, we may work only with functions which are smooth and such that the non-zero terms in | |7.1| l are finitely 
many. Then, since all of our results are quantitative and all the constants involved do not depend on i/i, by a density 
argument we may lower the regularity of -0 requiring only certain norms depending on the initial data of ^ to be 
finite. 

^Note that in the subextreme range T becomes spacelike in the region bounded by the two horizons, which 
however coincide in the extreme case. 
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We know by Proposition 5.2.1 that Jj^nf^ is non-negative definite. 



A.2 



However, we need to know 
for the definition of ^) is 



the exact way depends on d^. Note that ^„ (see Appendix 

strictly positive and uniformly bounded. However = —^g{T,T) ~ ^D, which clearly vanishes 
(to second order) at the horizon. Therefore, (A.l) implies 



where the constants in ^ depend on the mass M and Sg. Note that all these relations are invariant 
under the flow of T. 

On V.^, the energy estimate 



2) degenerates with respect to the transversal derivative dr^. It 



is exactly this that does not allow us to use estimate (8.2) to obtain the boundedness result for the 



waves in the whole region TZ. However, if we restrict our attention to the region where r > tq > M 



then commuting the wave equation with T and estimate (8.2) in conjunction with elliptic and 



Sobolev estimates give us the boundedness of in this region. This result is not satisfactory since 
it provides no information about the behaviour of waves on the horizon and so it is not sufficient 
for non-linear stability problems. 



9. MORAWETZ AND X ESTIMATES 

In view of the absence of the redshift along 7^+, any proof of uniform boundedness of solutions to 
the wave equation essentially relies on the dispersion of waves. Therefore, we must derive estimates 
which capture these dispersive properties. The first result we prove in this direction is that there 
exists a constant C which depends on M and Eg such that 



where the weight x — x(^) degenerates only at "H^, at r = 2M (photon sphere) and at infinity. 
The degeneracy at the photon sphere is expected in view of trapping. In particular, as we shall see, 
such an estimate degenerates only with respect to the derivatives tangential to the photon sphere. 
This degeneracy can be overcome at the expense of commuting with T. The degeneracy at will 
be overcome in Section [12] by imposing necessary conditions on the spherical decomposition of ip. 

9.1. The Vector Field X. We are looking for a vector field which gives rise to a current whose 
divergence is non-negative (upon integration on the spheres of symmetry) . We will be working with 
vector fields of the form X = f [r*) dr* (for the coordinate system {t, r*) see Section [2]). 

9.1.1. The Spacetime Term . We compute 

where /' = ^^^ip- and H — ^ '^^dr^ ■ Note that all the derivatives in this section will be considered 
with respect to r* unless otherwise stated. 

Unfortunately, the trapping obstruction does not allow us to obtain a positive definite current 
so easily. Indeed, if we assume that all these coefficients are positive then we have 



D 2D r \ D rj \D r ^ 

However, the quantity — changes sign exactly at the radius Q of the photon sphere. There- 
fore, there is no way to make all the above four coefficients positive. For simplicity, let us define 

P{r) {r~M){r~2M) 
= --n ■ r + U = . 

9.2. The Case I > 1. We first consider the case where ip is supported on the frequencies I > 1. 
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9.2.1. The Currents J^'^ and K^'^ . To assist in overcoming the obstacle of the photon sphere we 

shall introduce zeroth order terms in order to modify the coefficient of (StV')^ and create another 
which is more flexible. Let us consider the current 

jXMM,^ ^J^+h^ (r) VV^^ + h2 (r) {V ^w) , 

where /ii, /12, w are functions on Ai. Then we have 

^K^+ hi (v'^?AVaV')+ (v^/ii + 2/12^^) ^/-v^v + (v^/i2V^w + h2 (agw)) 

By taking h2 = = 2G,w — ~G we make the coefficient of ipV ^i) vanish. Therefore, let us 
define 

(9.1) J^'i = J^ + 2GV(V^V)-(V^G)^2^ 

Then 

Therefore, if we take G such that 



4 2r 



then 



r'^ V 4:D r D ■ r \ 2D ■ r 2r^ 



Note that, since / must be bounded and /' positive, — /"' must become negative and therefore has 
the wrong sign. Note also the factor D at the denominator which degenerates at the horizon. The 
way that turns out to work is by borrowing from the coefficient of {dr-'ip)'^ which is accomplished by 
introducing a third current. Note that this current exploits an algebraic property of the trapping 
in the extreme case (see the computation (9.4) of the term / below). 

9.2.2. The Current J^'^ and Estimates for K^''^. We define 

where /3 will be a function of r* to be defined below. Since Div (9^* ) = ^ + ^ , we have 

=(J(a,..^ + /3V)' + ^iy^l'- 



If" ( D\ f ( , , 2D D'\ f f D' D" 



Note that the coefficient of /i/;^ is independent of the choice of the function /3. Let us now take 
/3 = ^ — ^2^^i , where x = r* — a — ^/a and a > a sufficiently large number to be chosen 
appropriately. As we shall see, the reason for introducing this shifted coordinate x is that we want 
the origin a; = to be far away from the photon sphere. We clearly need to choose a function / that 
is stricly increasing and changes sign at the photon sphere. So, if we choose this function (which 
we call r ) such that (/")' = /« (r* = 0) = 0, then 

p.^ I in'" X in" in' ^ i x^~a^ 

AD a^+x^ D {a^ + x^f D 2D(x^ + a^f' 
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If we define 



X" = rdr> and / = ( ^ - j r then 




x-^ ''[^] - ^ + /sv-) + m\ +iF + i) r. 

9.2.3. Nonnegativity of . We have the following important proposition 

Proposition 9.2.1. There exists a constant C which depends only on M such that for all solutions 
ijj to the wave equation which are supported on the frequencies I > 1 we have 



We compute the term / 

Therefore, / > and vanishes (to second order) at the horizon and the photon sphere. Note now 
that the term F is positive whenever x is not in the interval [—a, a]. On the other hand, the 
photon sphere is far away from the region where x £ [—a, a], so / is positive there. However, / 
behaves like pi and so it is not sufficient to compensate the negativity of F. That is why we need 
to borrow from the coefficient of ^ip which behaves like ^. Indeed, the Poincare inequality gives us 
/§2 ^4''^ < /g2 lyV'l^ and, therefore, it suffices to prove that 2^-^ + F > for all a; G [—a, a] or, 
equivalently, r* G [y/a,2a + ^/a]. But 

^ 1 x^ — ^ x^ ~ a'^ 

~ D2{x^+a^f ~ "2(.t2 +a2)3 

with 1 < Ac 1 as a — >■ +oo. Moreover, since r* > r for big r, by taking a sufficiently big we 
have 

2^—^ > 26. ■' 



with 1 > (5q — 5- 1 as a — > +oo. Therefore, we need to establish that 

^ (q;2 - x^) {x + a + yS)^ da 

for all X € [—a, a]. In view of the asymptotic behaviour of the constants Aq,, Sa it suffices to prove 
that the right hand side of the above inequality is strictly less than 1. 

Lemma 9.2.1. Given the function /" (r*) = /" (r* (x)) defined above, we have the following: If 
-a<x<0 then /" > ftj;. Also, if 0<x <a then f" > - ^. 

Proof. For —a<x<Owe have 

J-a-^ a;2 + a2 x^ + 2a2 2a2 

Now for < X < a we have 

.X ^ n-x ^ nx ^ 

/" (a;) > / :^dx > / + / ^5 T^dx 

J -a X + " J -a X + J-x + 01 



^'^Note that the role of the parameter a as an index in our notation in this section is to emphasise that the 
corresponding tensors (i.e. functions, vector fields, etc.) depend on the choice of a. 
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□ 

Consequently, if —a < x < and x = —Xa then < < 1 and 

ia-x){x + a + q2 _ (1 - A) (l + A + g-s) ^ ^ (1 - A) (1 + A)^ 2 _9_ 

2(a:2 + a2f 2 (A2 + if 2 (A2 + 1)' 3 10' 

since c?q — !■ 1 as a — > +oo. Similarly, if < x < a and x — \a then < A < 1 and 

^ a (a2 _ a;2j (^^ + a + a {a^ ~ x"^) {x + a + ^/af _ ~ (l - A2) (1 + A)^ 9 

2 {x^ + a^f {a^ + 2ax - x-^) ^ 2(a;2 + a2)^ ~ " 2 (1 + A2)^ '^10' 

since da 1 as a +oo. □ 

Note, however, that although the coefficient of ip does not degenerate on the photon sphere, the 
coefficient of the angular derivatives vanishes at the photon sphere to second order. Having this 
estimate for -0, we obtain estimates for its derivatives (in Section 9.7 we derive a similar estimate 
for ■0 for the case 1 = 0). 

Proposition 9.2.2. There exists a positive constant C which depends only on M such that for all 
solutions Tp of the wave equation which are supported on the frequencies I > 1 we have 



Proof. From (9.3 1 we have that the coefficient of ip does not degenerate at the photon sphere. The 
weights at and infinity are given by the Poincare inequality. Commuting the wave equation 
with T completes the proof of the proposition. □ 

9.2.4. The Lagrangian Current Lj^. In order to retrieve the remaining derivatives we consider the 
Lagrangiarp] current 

where / = ^^Di. 

Proposition 9.2.3. There exists a positive constant C which depends only on M such that for all 
solutions tjj of the wave equation which are supported on the frequencies I > 1 we have 



IWl'j (^Div(L^) + C^if^°'2 [TV]j. 



Proof. We have 

Div(L;^) =/V^VV^V + V^/V-V^V; 



/I? _ L>3/2 , \/Z?_ ,., ID,., D 



vtKj ryt'J rft\.} ^ fyt'-t rft'-t 



where e > is such that t\/lD < |. Therefore, in view of Proposition 9.2.2 we have the required 
result. □ 

Proposition 9.2.4. There exists a positive constant C which depends only on M such that for all 
solutions ijj of the wave equation which are supported on the frequencies I > 1 we have 

I (^*^)^ + ll (^-V.)^ + ^ \nf + ^^^) < I (Div(X/) + C±K-^^^ [TV]) . 

Proof. Immediate from Propositions |9.2.T] and [9?273l □ 

^^The name Lagrangian comes from the fact that if ip satisfies the wave equation and C denotes the Lagrangian 
that corresponds to the wave equation, then C{tp, dip, g~^) = g'^'^ dfitpd^tp = Div(</)V^i/') 
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9.2.5. The Current . In case we allow some degeneracy at the photon sphere we can obtain 

similar estimates without commuting the wave equation with T. This will be very useful for 
estimating error terms in spacetime regions which do not contain the photon sphere. We define the 
function such that 



1 



(r*)2 + r 



f{r* = 0) = 0. 



Proposition 9.2.5. There exists a positive constant C which depends only on M such that for all 
solutions ijj of the wave equation which are supported on the frequencies I > 1 we have 



1 

C 



§2 



< 



§2 



CK 



where X = f 8^* and the current J^' is as defined in Section 



9.2.1 



Proof. Note that the coefficient of ilP in vanishes to first order on (see also Lemma 9.4.1 ) 

and behaves like ^ for large r. Note also that the coefficient of {dr'ip)^ converges to AP (see again 
Lemma 9.4.1 ). The result now follows from Proposition 9.2. 1| □ 

In order to retrieve the 9t-derivative we introduce the current L^^f = /I'^-^V^V; where h'^ is such 
that: 



1 



(r*)2 + l 



for M <r <ro < 2M, ft,'' < for < r < 2M, 



h'^ for r = 2M to second order, ft'' < for 2M < r < ri, ft'' = - ^ for ri < r. 



Proposition 9.2.6. There exists a positive constant C which depends only on M such that for all 
solutions tp of the wave equation which are supported on the frequencies I > 1 we have 



1 
C 



§2 



P • (r - 2M) 



< 



§2 



Proof. We have as before 

Div(<) = -^(9tV)' + ^{dr-i^f + h" IWI^ + (^V'(9..V)- 

Since the coefficient of ^p{dr''ijj) vanishes to first on the horizon, the Cauchy-Schwarz inequality and 
Proposition |9 . 2 . 5] imply the result. □ 

Finally, we obtain: 

Proposition 9.2.7. There exists a positive constant C which depends only on M such that for all 
solutions Tp of the wave equation which are supported on the frequencies I > I we have 

1 f /P- (r-2M),^ 1,^ ,^2 P ■ (r -2M) ,^ ^,2 (r-Af),, 

U Jg2 \ r° 

< 



/ (Div(L^') + if^'-i[^]+Cif^°''[V'] 

JS2 ^ 



Proof. Immediate from Propositions |9.2.T] and [93761 



□ 
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9.3. The Case / = 0. In case I = the wave is not trapped. Indeed, if we consider the vector field 
X" = fdr- with f then we have: 

Proposition 9.3.1. For all spherically symmetric solutions %[) of the wave equation we have 



Proof. Immediate from the expression of and the above choice of / = /°. 



□ 



9.4. The Boundary Terms. The positivity of the bulk terms shown so far will be useful only if 
we can control the arising boundary terms. 

9.4.1. Estimates for J^n^. 

Proposition 9.4.1. Let X = fdr* where f = f{r*) is hounded and S he either a 5*0(3) invariant 
spacelike (that may cross ) or a SO(S) invariant null hypersurface. Then there exists a uniform 
constant C that depends on M , S and the function f such that for all we have 



< c 



Proof. We work using the coordinate system (w, r, 6*, (/>). First note that X — fdr* — fd^ + f ■ Ddr 
Now 



= ifn-) {d^ijjf + Dfn" (a.V) + D 



lofn" - Ifn' - lofn" + fn"^ 
2 2 2 



-Dfn" - -Dfn'' - -fn' 
2 2 2 



The result now follows from the boundedness of /. 
9.4.2. Estimates for J^°'* ng, i ^ 1,2. 

Proposition 9.4.2. There exists a uniform constant C that depends on M and S such that 



□ 



J. 



5 



where S is as in Proposition 9.4-1 
Proof. It suffices to prove 



2G"^(V^7A)-(V^G")^/.2 + 



in' 

D 



Js 



For we first prove the following lemma that is true only in the case of extreme Reissner-Nordstrom 
(and not in the subextreme range). 



Lemma 9.4.1. The function 



F 



1 



1 



Z?((r*)2 + 1) 
is bounded inTZU li^ . 

Proof. For the tortoise coordinate r* we have r* (r) = r + 2M ln(r - M) - + C. Clearly F ^ 
as r — > +0O. Moreover, in a neighbourhood of T-L^ 

.2 1 



F 



{r - M)2 
This implies the required result. 



_M^+^2+4M2(ln(.-Af))2 



< 00. 



□ 
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An immediate coUorary of this lemma is that the fmictions Gf = and ''^j^ are bounded 
in 7^ U Furthermore, for M < r we have (/")' - ^, D - 1 and drD - ^. Therefore, 

= r^VrG"" is also bounded. Finally, WyG" = 0. The above bounds and the first Hardy 
inequality complete the proof of the proposition. □ 

9.5. A Degenerate X Estimate. We first obtain an estimate which does not lose derivatives but 
degenerates at the photon sphere. 

Theorem 9.1. There exists a constant C which depends on M and Eq such that for all solutions 
ijj of the wave equation we have 



(9.5) 



K(0/ 



#^2 



Proof. We first decompose '0 as 

We apply Stokes' theorem for the current 



LI [0; 



in the spacetime region 7?,(0, r) and use Propositions 9.2.7|and 9.4.2 We also apply Stokes' theorem 
in TZ{Q,t) for the current [-00] and use Proposition 9.4.1 and by adding these two estimates we 
obtain the required result. 

□ 



Remark 9.1. One can further improve the weights at infinity. Indeed, if we apply Stokes' theorem 
in the region {r > R} for sufficiently large R for the current 



where /i = 1 — /2 = 2^7^ and 6 > 0, then an easy calculation shows that there exists a 
constant Cs > such that for r > R 



'^''J^>Cs(r-'-'{dttPf 



„-3-<5 



0^) 



The above remark and theorem and the identity d^* = d^, + Ddr imply statement (2) of Theorem 
[l]of Section |4) 

9.6. A non-Degenerate X Estimate. We derive an estimate which does not degenerate at 
the photon sphere but requires higher regularity for -0. 

Theorem 9.2. There exists a constant C which depends on M and Sq such that for all solutions 
ip of the wave equation we have 



(9.6) 



K(0,r)\ ^ J •'^0 



Proof. We again decompose ip as tp = 0>i + ipQ and apply Stokes' theorem for the current 



in the spacetime region TZ{0, t), where G is the constant of Proposition 9.2.4 and use Propositions 
and 



9.2.4 



9.4.2 



Finally, we apply Stokes' theorem in TZ{0, r) for the current Jf^ [0o] and use 



Proposition |9.4. l] Adding these two estimates completes the proof. □ 

The proof of statement (3) of Theorem [l] of Section [4] is immediate from the above proposition. 
Theorem 19. II and Remark 1 9. II 
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Remark 9.2. All the above estimates hold ifS^ is replaced by E^.. Indeed, the results of Section 9. 
allow us to bound the corresponding boundary terms in this case. Moreover, the additional boundary 
term on X"*" can be bounded by the conserved T-flux. 

9.7. Zeroth Order Morawetz Estimate for tp. We now prove weighted estimates of the 



wave ip itself. In Section 9.2.3 we obtained such an estimate for tpyi. Next we derive a similar 
estimate for the zeroth spherical harmonic -00- We first prove the following lemma 

Lemma 9.7.1. Fix R > 2M . There exists a constant C which depends on M , Sq and R such that 
for all spherically symmetric solutions ip of the wave equation 

J{r=i?}nK(0,T) JSo 



Proof. Consider the region J-{Q, t) = {7^(0, t) n {M < r < R}} . By applying the vector field X 
dr" as a multiplier in the region J-(0, r) we obtain 



where Ujr denotes the unit normal vector to {r = R} pointing in the interior of J-. 



















Q 




II 





1+ 



In view of the spatial compactness of the region T , the corresponding spacetime integral can 
be estimated using Propositions 



estimated using Proposition 



9.4.1 



9.3.1| and |9.4.1[ The boundary integrals over Eg and E,- can be 
which completes the proof. □ 



Moreover, for spherically symmetric waves V' we have [V'I'^^f 



Consider now the region Q = 'R,V\{R<r^ . 

Proposition 9.7.1. Fix R > 2Af . There exists a constant C which depends on M , Eg and R such 
that for all spherically symmetric solutions tjj of the wave equation 

1 1.2 



r<CR 



{r=R} 



Proof. By applying Stokes' theorem for the current J'^'^[^] in the region Q{Q,t), where again 
X = dr' (and, therefore, / = 1), we obtain 



IQ{0,t) J{r=R}ng(0,T) JSo 

where we have used Proposition |9.4.2] to estimate the boundary integral over E,- n Q. Note that 
again Ug denotes the unit normal vector to {r = R\ pointing in the interior of Q. For / = 1 we 
have K ^''^[%li\ = I^j^, where 7 > and / ;;ir for r > i? > 2M. If G is the function defined in 
Section 



9.2.1 



then G 



1^ and therefore for sufhciently large R we have dr»G < 0. Then, 



^9,*, by applying Cauchy-Schwarz for the second term on the right hand side (and 



Since ng 

since the third term is positive) we obtain 



e 
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for a sufficiently small e. Lemma [9. 7. 1 1 completes the proof. □ 

It remains to obtain a (weighted) LP' estimate for ijj in the region J- . 

Proposition 9.7.2. Fix sufficiently large R > 2M . Then, there exists a constant C which depends 
on M , So o.'fT'd R such that for all spherically symmetric solutions of the wave equation 



Proof. We apply Stokes' theorem for the current 

where H — {r ~ Af)^. All the boundary integrals can be estimated using Propositions 9.4.2 and 
[9jn Note also that 

Since OgH ~ D for M < r < R and j; is bounded, the result follows in view of the spatial 



compactness of and Proposition 9.3.1 □ 



We finally have the following zeroth order Morawetz estimate: 

Theorem 9.3. There exists a constant C that depends on M and Eq such that for all solutions 
of the wave equation we have 

(9-7) / ^i^'<C f Jj[#.5,„ 



n 



Proof. Write = ipa + V-'>i a^nd use Propositions 9.2.1 9.7.2 and 9.7.1 



□ 



Clearly, the lemma used for Proposition 9.7.1 



could have argued by averaging the X estimate (9.51 in R, which would imply that there exists a 



holds stricly for the case 1 = 0. In general, one 



value Ro of r such that all the derivatives are controlled on the hypersurface of constant radius Rq. 
This makes our argument work for all -0 without recourse to the spherical decomposition. 



Remark 9.3. The key property used in Proposition 9.7.1 
negative, something not true for larger powers of - 



is that the d' Alembertian of - is always 



Note that this unstable behaviour of ^ is 
expected since it is the static solution of the wave equation in Minkowski. 



9.8. Discussion. The current of Section 9.2.2 was first introduced in pT] and subsequently in [23], 
where a non-degenerate X estimate is established for Schwarzschild. In [21], the authors show 
that for each fixed spherical number / there is a corresponding "effective photon sphere" centred at 
r — — 7;. It is also shown that lim/_j._|_co —7/ — >■ 3A/. In our case, our computation ( |9.4[ ) shows that 
all "effective photon spheres" coincide with the photon sphere, i.e. for all I we have —7/ — 2M. 
Note also that in [23], one needs to commute with the generators of the Lie algebra so(3). Moreover, 
the boundary terms could not be controlled by the flux of T but one needed a small portion of the 
redshift estimate. In our case the above geometric symmetry of the trapping and the asymptotic 
behavior of r* allowed us to completely decouple the dispersion from the redshift effect. For a nice 
exposition of previous work on X estimates see [2^ . 



10. The Vector Field N 

It is clear that in order to obtain an estimate for the non-degenerate energy of a local observer we 
need to use timelike multipliers at the horizon. Then uniform boundedness of energy would follow 
provided we can control the spacetime terms which arise. For a suitable class of non-degenerate 
black hole spacetimes, not only have the bulk terms the right sign close to but they in fact 
control the non-degenerate energy. Indeed, in [3S] the following is proved 
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Proposition 10.0.1. Let be a Killing horizon with positive surface gravity and let V be the 
Killing vector field tangent to H^. Then there exist constants B > and a timelike vector field 
N which is (j)^ -invariant (i.e. CyN = 0) such that for all functions "0 we have 



(10.1) 
onn+. 



The construction of the above vector field does not require the global existence of a causal Killing 
vector field and the positivity of the surface gravity suffices. Under suitable circumstances, one can 
prove that the N flux is uniformly bounded without understanding the structure of the trapping 
(i.e. no X or Morawetz estimate is required). 

However, in our case, in view of the lack of redshift along "H^ the situation is completely different. 
Indeed, we will see that in extreme Reissner-Nordstrom a vector field satisfying the properties of 
does not exist. Not only will we show that there is no 0!^-invariant vector field 
on T-L^ but we will in fact prove that there is no 0^-invariant tinielike vector 



Proposition 10.0.1 



N satisfying (10.1 



field N such that 



if'' [i^] > 



on 7^"*". Our resolution to this problem uses an appropriate modification of and Hardy inequal- 
ities and thus is still robust. 



10.1. The Effect of Vanishing Redshift on Linear Waves. Let N = N'" (r) + N'' (r) dr be 

a future directed timclikc (7i)!r-invariant vector field. Then 



where the coefficients are given by 
(10.2) 

F,, = (5,iV") ,Frr = D 



2 7 



N^'D' 1 2N'^ 

,Fy = --{drNn,F,r^D {drN^ - — 



where D' 



^ . Note that since g {N, N) 



27V" iV, g{N,T) 



N^, and so 



{r — M) can not be zero (otherwise the vector field N would not be timelike). Therefore, looking 



back at the list ( 10.2 ) we see that the coefficient of {dr^) vanishes on the horizon 7^+ whereas the 



coefficient of d^tpdr^ is equal to which is not zero. Therefore, K [ip] is linear with respect 

to dri^ on the horizon "H^ and thus it necessarily fails to be non-negative definite. This linearity 
is a characteristic feature of the geometry of the event horizon of extreme Reissner-Nordstrom and 
degenerate black hole spacetimes more generally. This proves that a vector field satisfying the 
properties of Proposition 1 1 . . l] does not exist. 



10.2. A Locally Non-Negative Spacetime Current. In view of the above discussion, we need 
to modify the bulk term by introducing new terms that will counteract the presence of dy-ipdri^. 
We define 



J = J'^''' - J 



(10.3) 

where /i is a function on A4. Then we have 



N 



h (r) VV^V 



K = K 



N.h 



- K'' + (V^h) ipv^^j + h (V^^VaV) , 



provided is a solution of the wave equation. Let us suppose that h (r) ~ 



_ N''(M) 



M 



Then 



(10.4) 



K^'''[^] = K^i^p] + h (V^T/^Va^/-) = if^N + h (2dy^dr^ + D {dri'f + 



(a„V)' + [Frr + hD] {dr^f + 



Ifijf + [F„r + 2h] {d,4'dri^) . 
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Note that by taking h to be constant we managed to have no zeroth order terms in the current K. 
Let us denote the above coefficients of [da'4'dh^) by Gab where a, 6 e {v, r, ^} and define the vector 
field N in the region M < r < ^ to be such that 

(10.5) N'^ir) = 16r, iV'^(r) = -^r + M 

and, therefore, /i = — |. Clearly, N is timclikc future directed vector field. We have the following 



Proposition 10.2.1. For all functions V', the current K-'^'~2\^p] defined by (10.4) is non-negative 
definite in the region An = {M < r < and, in particular, there is a positive constant C that 

depends only on M such that 



K 



Proof. We first observe that the coefficient Gy of dyipdrip is equal to Gy — 1. Clearly, Gw = 16 
and Grr is non-negative since the factor of the dominant term D' (which vanishes to first order on 
is positive. As regards the coefficient of the mixed term we have Gyr — 16-D + = ei • £2, 
where ei = -s/I?, £2 — ^QVD + 2. We will show that in region An we have ef < Grr: ^2 < G^y. 
Indeed, if we set A = 4f we have 



4 < Grr ^ (1 - A) < (1 - A) - Aj - A (^-- + Aj ^ A < -, 

which holds. Note also that ef vanishes to second order on whereas Grr vanishes to first order. 
Therefore, in region we have Grr ~ t\ ^ \/D. Similarly, 



'2 < G„„ ^ 16 









+ 2 







7 

< 16 4^ A > 



which again holds. The proposition follows from the spatial compactness of An and that -\- ah-\- 
6^ > for all a, 6 e M. □ 

10.3. The Cut-off i5 and the Current J^'^' ^. Clearly, the current K^'~^ will not be non- 
negative far away from "H^ and thus is not useful there. For this reason we extend N'",N^ such 
that 

N" (r) > for all r > M and N" (r) 1 for all r > 

N'' (r) < for all r > M and N'' (r) = for all r > 

N as defined is a future directed timelike (/)^-invariant vector field. Similarly, the modification term 

n 

in the current ' ^ is not useful for consideration far away from 7^+. We thus introduce a smooth 
cut-off function S : [A/,-foo) ^ R such that S{r) = l,r e [M, ^M] and S{r) =0,r e [^,+00) 
and consider the currents 

(10.6) jH^''-^- = j; - ^<5^V,^, K^'''-'^ = 

We now consider the three regions An,BntCn 
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In region Cn ~ {r > ^} where (5 = and N = T, we have K^'^'~i — 0. However, this 
spacetime current, which depends on the 1-jet of '0, wih generally be negative in region Bn and 
thus will be controlled by the X and Morawetz estimates (which, of course, are non-degenerate in 
Bn). 

We next control the Sobolev norm IIV'll^i(s^) ^ (dvi^)'^ + {drip)'^ + I^V'I^- Since N and 
are timelike everywhere in TZ and since ^Ar,^„ are positive and uniformly bounded, (A.l) of 



Appendix A. 2 implies Jf^ [ip]n'^ ^ (dvip) + (drtp) + lyV-"! i where, in view of the 0!^-invariance 
of and iV , the constants in ~ depend only on M and Eq- Therefore, it suffices to estimate the 
flux of N through St-. We have: 

Proposition 10.3.1. There exists a constant C > which depends on M and Sq such that for all 
functions ip 

(10.7) / J^mn^^<2[ j;!''-^^]n>^ + C [ JJ[VK. 



Proof. We have 

> J^n^ - e5{d,i^f - e5{dr^f - -^^^ > ^ J^n^ - -^^^ 

for a sufficiently small e. The result follows from the first Hardy inequality. □ 

Corollary 10.1. There exists a constant C > which depends on M and Sq such that for all 
functions ip 

(10.8) / Jl!'''~^i,]n^ <C [ J^[^]n^. 

Proof Note that 

and use the first Hardy inequality. □ 

10.4. Lower Estimate for an Integral over T-L^ . As regards the integral over "H+j we have 
Proposition 10.4.1. For all functions tp and any e > we have 



Proof. Recall our convention = T . On "H^ we have (5 = 1. Therefore, J^]'"' ^ = ^7!^+ ^ 
^ipdyip. However, 

From the first and second Hardy inequality we have 

n+n's JT, JY, JE JE 

which completes the proof. □ 



W+ JW+ JT.^ JS^ 

where Ce depends on M , Sq and e. 
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11. Uniform Boundedness of Local Observer's Energy 
We have ah tools in place in order to prove the following theorem 

Theorem 11.1. There exists a constant C > which depends on M and Sq such that for all 
solutions tp of the wave equation 



(11.1) 



Proof. Stokes' theorem for the current Jf 



N,S,- 



in region TZ{0,t) gives us 



K 



n 



J, 



N,5.-i 



H + 



J, 



First observe that the right hand side is controlled by the right hand side of (10.8). As regards the 



left hand side, the boundary integrals can be estimated using Propositions 10.3.1| and 10.4.1 The 
spacetime term is non- negative (and thus has the right sign) in region An, vanishes in region Cn and 
can be estimated in the spatially compact region Bn (which does not contain the photon sphere) 
by the X estimate (9.5) and Morawetz estimate (9.7). The result follows from the boundedness of 
T-flux through S^. □ 

Corollary 11.1. There exists a constant C > which depends on M and Sq such that for all 
solutions ip of the wave equation 



(11.2) 



K 



Theorem [2] of Sectio n [4| is then implied by ( |11.1[ ), Proposition 10.2.1 and the above corollary. 



Note also that Corollary 11.1 and estimate ( |9.6p give us a spacetime integral where the only weight 
that locally degenerates (to first order) is that of the derivative tranversal to . Recall that in 
the subextreme Reissner-Nordstrom case there is no such degeneration. In the next section, this 



degeneracy is eliminated provided ■00 = 0. This condition is necessary as is shown in Section 12.1 



One application of the above theorem is the following Morawetz estimate which does not degen- 
erate at V.^. 

Proposition 11.0.2. There exists a constant C > which depends on M and Eg such that for all 
solutions Tp of the wave equation 



(11.3) 



Proof. The third Hardy inequality gives us 



D{[d,^f + {drij)\] 



where C is a uniform positive constant that depends only on M . The integral over Bn of ip'^ can 

and the last integral on the right hand side can be estimated using (9.5) 
and Corollary 1 11.1 1 



be estimated using (9.7 



and Propositions 1 1 . 2TT 



□ 



The above proposition in conjunction with Remark 9.2 and Theorem 9.3 implies statement (1) 
of Theorem H] of Section [H 

Remark 11.1. Note that all the above estimates hold if we replace the foliation with the foliation 
Y.T- which terminates at 1+ since the only difference is a boundary integral over Z+ of the right sign 
that arises every time we apply Stokes' theorem. The remaining local estimates are exactly the 
same. 
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12. Commuting with a Vector Field Transversal to 



As we have noticed, the estimate ( 11.2 1 degenerates with respect to the transversal derivative dr 
to T-L^ . In order to remove this degeneracy we commute the wave equation □gV' — with dr aiming 
at additionaUy controUing all the second derivatives of -0 (on the spacelike hypersurfaces and the 
spacetime region up to and including the horizon 7^+). Such commutations first appeared in [M] . 

12.1. The Spherically Symmetric Case. Let us first consider spherically symmetric waves. We 
have the following 

Proposition 12.1.1. For all spherically symmetric solutions ijj to the wave equation the quantity 

(12.1) ^-^+]^^ 

is conserved along . Therefore, for generic initial data this quantity does not decay. 

Proof. Since i/j solves Ogip — and since /^ip = we have dydr'4' + jjdv4' = and, since is 
tangential to Ti^ , this implies that drip +10^^ remains constant along Ti^ and clearly for generic 
initial data it is not equal to zero. □ 

By projecting on the zeroth spherical harmonic we deduce that for a generic wave ip either ip'^ 
or {dripY does not decay along Ti.^ completing the proof of Theorem [s] of Section |4j The nature 
of the above proposition (which turns out to be of fundamental importance for extreme black hole 
spacetimes) will be investigated in great detail in the companion papers [21 H] , where, in particular a 
similar non-decay result is shown to hold even for "generic" waves whose zeroth spherical harmonic 
vanishes. 

Clearly, this proposition indicates that if one is to commute with dr then one must exclude the 
spherically symmetric waves and thus consider only angular frequencies I > 1. Indeed, we next 
obtain the sharpest possible result. 

12.2. Commutation with the Vector Field 9,.. We compute the commutator [□g,c}r]- First 
note that dr] ip = f ^V'- Therefore, if i? = + = ^ then 

(12.2) [Dg, dr] ^ = -D'drdr^ + \d,ip - R'dr^ + -^^A, 



r 



where i?' = ^ . 

12.3. The Multiplier L and the Energy Identity. For any solution tjj of the wave equation we 
have complete control of the second order derivatives of ip away from "H^ since 

(12.3) \\dadt^fm^^n{M<n,<r}) <G j Jlm<, +C f jj [T^] 

where C depends on M, ro and Eq and a,b £ {v,r,y/}. Note that we have used elliptic estimates 
since T is timelike away from "H"*" (the zeroth order terms can be estimated by the first Hardy 
inequality). Similarly, away from , we have control of the bulk integrals of the second order 



derivatives since (9.7) and local elliptic estimates imply 



(12.4) \\dad,ij]\l.^^^o.r)n{M<ro<r<r,<2M}) <C jj [V'] + C jj [T^] • 

J So -'So 

The above estimates will be required for bounding several spacetime error terms away from T-i^ . 

In order to estimate the second order derivatives of "0 in a neighbourhood of T-L^ we will construct 
an appropriate future directed timelike 0!^-invariant vector field 

L = L''d^+L''dr, 

which will be used as our multiplier. Since we are interested in the region M < r < ro < ri, L 
will be such that _L = in r > ri and L timelike in the region M < r < ri. Note that ro,ri are 
constants to be determined later on. The regions A, B are depicted below 
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For simplicity we will write TZ instead of TZ (0, r), A instead of A (0, t), etc. The "energy" identity 



for the current J^' [drip] is 



(12.5) 



's^ Jn Jh+ JSo 

The right hand side is controlled by the initial data and thus bounded. Also since L is timelike in 



the compact region A we have from (|A.1|) of Appendix |A.2 
(12.6) 



where the constants in ^ depend on M, Sq and L. Furthermore, on we have 



(12.7) 



'«+ - ,\^ v^; -r , 2 

Therefore, the term that remains to be understood is the bulk integral. Since drip does not satisfy 
the wave equation, we have 

[drip] ^ [dr^j] + [dr^] = [Or^j] + (Dg {dr1p)) L {drip) . 

We know that 

[drip] - F„„ {d^drlpf + Frr {drdrlpf + [fdr1p[^ + F^r (d^Orlp) (drdrlp) 



where the coefficients are as in (10.2 1. In view of equation (12.2) we have 

£^ [drip] = - D'i'' {drdrlpf - DT (drdrlp) (drdrlp) " i^'L" (dM) (drlp) 



(dydrip) (dyip) 



2—{dydrip)4\pJ- 



(drdrlp) (dyip) - R'U [drdrlp) {drip) 



2 — {drdrlP)H. 

r 



Therefore, we can write 

[dri^] =i/l {dydriP)'' + (drdrl/j)^ + H3 [fdrlP[^ + H4 (dydrl/j) (d,^) + H, (9,5,^) (5.V) 
+ He (drdrlp) (dyPj) + Hj{dydr1p)/t^1p + Hs{drdr1p)4Mp + Hg [drdrlp) [drdrlp) 
+ HiQ [drdrlp) [drip) , 

where the coefficients Hi,i — 1, 10 are given by 

^[drL"-) U 



(12. 



i?i = [drL-) , H2=D 
Hi 



r 



n' 1 

-^L\ H, = --[drU)., 



-2 — , H^ 



rR, He 



Hg = D[drU') - D'V - 2 



HiQ = — R' . 



2 — , i^s 
r 



Since L is a future directed timelike vector field we have L" [r) > and L*" (r) < near 'H'^. By 
taking drL^ [M) sufficiently large we can make Hi positive close to the horizon H+. Also since 
the term D vanishes on the horizon to second order and the terms R, D' to first order and since 
[M) < 0, the coefficient H2 is positive close to "H^ (and vanishes to first order on it). For the 
same reason we have HgD < ^ and (HgR) < ^ close to . Moreover, by taking ~drL^ [M) 
sufficiently large we can also make the coefhcient H^ positive close to such that Hg < j^. 
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Indeed, it suffices to consider U such that — ^ j*^-* < — ^''"^25^^'' ^'^'^ then by continuity we have 
the previous inequahty close to "H^ . Therefore, we consider M < tq < 2M such that in region 
A = {M < r < To} we have 

L" > 1, drL" > 1, ~r > 1, i/i > 1, H2 > 0, i?3 > 1, 
(12.9) H'i H9 9 H9 H-^ 

8 10 ' - 10 ' ^ ^ - 10 ' 10 

Clearly, ro depends only on M and the precise choice for L close to H'^. In order to define L 
globally, we just extend and such that 

L" > for all r <ri and = for all r >ri, 
— L*" > for all r < ri and — for all r > ri, 

for some ri such that tq < ri < 2Af . Again, ri depends only on M (and the precise choice for L). 
Clearly, L depends only on M and thus all the functions that involve the components of L depend 
only on M. 

12.4. Estimates of the Spacetime Integrals. It suffices to estimate the remaining 7 integrals 
with coefficients Hi's with i = A, 10. Note that all these coefficients do not vanish on the horizon. 
We will prove that each of these integrals can by estimated by the N flux for -0 and Tij) and a small 
(epsilon) portion of the good terms in K^ldrii]- First we prove the following propositions. 

Proposition 12.4.1. For all solutions tJj of the wave equation and any positive number e we have 
where the constant depends on M , Sq and e. 

Proof. By applying the third Hardy inequality for the regions A, B we obtain 

[dri^f <C f {dr^jf + C f D [{d^dr^Pf + {Ordri^f] , 
A JB JAUB 

where the constant C depends on M and Eq- Moreover, since CD vanishes to second order at H^, 
there exists with M < < tq such that in the region {M < r < r^} we have CD < eHi and 
CD < eH2. Therefore, 

[dri^f <C[ {dr^f + I eHiid^dr^f + ei?2 (^rC*. V)' +/ D [{dydri^f + (9,9,0)'] • 

Jb J{M<r<r,} J{r^<r<ri} 



The first integral on the right hand side is estimated using (9.5) and the last integral using the local 
elliptic estimate (12.4) since T is timelike in region {r^ < r < ri} since M < r^. The result follows 
from the inclusion {M < r < r^} C A and the non-negativity of Hi, i = 1, 2 in ^. □ 



Proposition 12.4.2. For all solutions ip of the wave equation and any positive number e we have 

<cJ + ^ J^^M< 

J'So 

where the positive constant Cg depends on M , Eg and e. 
Proof. Integrating by parts gives us 



u+ 



SqUSt 



(9,0)(a,V) = - / 0(5,9,0)+/ 0(9,V)-/ 0'(9,0). 
Since solves the wave equation, it satisfies —9,9,0 = ^^(9,0) + |<j2i0 on Ti^. Therefore 
0(9„9,0) =j- [ 0(9,0) + 11 0(^0) 
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All the integrals on the right hand side can be estimated by J^^ J^H']'''^'^^ using the second Hardy 

inequality and Furthermore, J-^+f^^ 4'{dr4') < In+nT, ^ In+ns (^'•V')^- From the first and 

second Hardy inequality we have 



where C depends on M and Sq. In addition, from the second Hardy inequality we have that for 
any positive number e there exists a constant which depends on M, Eq ^nd e such that 



«+ns 



□ 



Proposition 12.4.3. For all solutions ijj of the wave equation and any positive number e we have 



H + 



where the positive constant Cg depends on M , Eg (^^d, e. 
Proof. In view of the wave equation on "H^ we have 



H+ 



H+ 



The first integral on the right hand side can be estimated using Proposition |12.4.2 For the second 
integral we have 



H + 



H + nT.0 



where, as above, e is any positive number and Ce depends on M, Eq and e. 



'St' 



□ 



Estimate for / H4 {dvdri^) (dvtp) 
Jn 



For any e > we have 



A 



H4 (dydrip) {dyi(;) 



'So' 



where the constant depends only on M,Eo and e. 



Estimate for / H^^ (dydrip) {dri^) 
Jn 



As above, for any e > 



[dydril)) [drip) 



A 



where m = max_4 e ^H^. Then from Proposition 12.4.1 (where e is replaced with ^) we obtain 

H5 {dydr^p) {dr^) 



(12.10) 

where Ce depends on M, Eq and e. 



A 



< / e{dydr^y+ / eHiid^dr^y + eH2{drdr^y 



Estimate for / Hq (drdri^) {dyx/j) 
Jn 
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Since Divdr = Stokes' theorem yields 



n 



TZ 



n 



He (dyip) (drip) dr ■ 



^6 {dv-ip) (drip) dr ■ n^. 



He {dv-ip) {dr'ip)dr ■ n^+. 



The boundary term over can be est imated by the A^-flux whereas the boundary integral over 
"H^ can be estimated using Proposition 12.4.2 If Q = drHe + -He, then it remains to estimate 



n 



Q (d^ij) (dr^) + He{dM){dri^)- 



n 



The first integral can be estimated using Cauchy-Schwarz, (11.21 and Proposition 12.4.1 The 
second integral is estimated by (12.10) (where He, is replaced with He). 

Estimate for / H-j (d^drtp) /^V' 
Jn 

Stokes' theorem in region TZ gives us 

Hr{dM)4^^ + I Hj{dr^j){4^d.,,p) Hj (dr^) {^ij)d, ■ 



n 



Hj (drip) dy ■ ns^. 



For the boundary integrals we have the estimate 



where depends on M, Sq and e. As regards the second bulk integral, after applying Stokes' 
theorem on §^(r) we obtain 



/ HrfdriP ■ fdyij < e f \fdri^f 

J A J A 



A 



where depends on M, Sq and e. Note that the second integral on the rig ht ha nd side can be 
bounded by J^^ J^[T%l:]n'^^. Indeed, we commute with T and use Proposition |ll.l| (Another way 
without having to commute with T is by solving with respect to /f^ip in the wave equation. As we 
shall see, this will be crucial in obtaining higher order estimates without losing derivatives.) 



We have 
Jn 



Estimate for / (drdrip) 4^ip 



n 



dr {HsH) (a,V) + / Hs im (^rV-) 



n 



n 



■H + 



The integral over is estimated in Proposition 12.4.3 Furthermore, the Cauchy-Schwarz in- 
equality implies 



Hs (f V) (drij) dr ■ ns. 



^-tHA 



Hs (yV ■ fdr^) dr ■ ns. 



s^nA 



where e is any positive number and Cg depends only on M, Sq and e. For the remaining two 
spacetime integrals we have 

2 



dr {H^m {dr^) 



-Hs (f V) 



-Hs lys.vr 



-drHsif^-fdri'). 



The conditions (12.91 assure us that \Hs\ < and thus the first integral above can be estimated. 
For the second integral we apply the Cauchy-Schwarz inequality. 
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Estimate for / Hg {dydrtp) (drdrip) 
Jn 

In view of the wave equation we have dydrip ~ —^drdrip— ^dy^ ~ ^drip — ^4^ip- Therefore, 



Hg (dydrlp) (drdrlp) = - Hg 



2 2 r 2 



Note that in A we have HgD < ^ and thus the first term on the right hand side poses no problem. 
Similarly, in A we have 

According to what we have proved above, the integral — ^ (dyip) (drdri^) can be estimated 
(provided we replace iJg with — ^). Similarly, we have seen that the integral — ^{drdril>)4^ip 
can be estimated provided we have Hg < which holds by the definition of L. 



We have 



Estimate for / Hiq (drdrip) {drip) 



Hio {drdri^) (dri^) 



n 



n 



n 



drH 



10 



-H 



10 



= / HiQ {drijj) dr ■ - 



Hio (drlpf dr • - / Hio {drtpf. 

'So J^r J'H + 

The second spacetime integral can be estimated using Proposition |12.4.1| whereas the boundary 
integral over E,- can be estimated using the A''-flux. Finally we need to estimate the integral over 

Ti.^ . It turns out that this integral is the most problematic. Since R' = D" + ^ we have 

R' (M) = ^2 and thus 

(12.11) Hxo{M) ^ -L"" {M)R! {M) ^ -U [M) ^ 



A/2 



> 0. 



In order to estimate the integral along H"*" we apply the Poincare inequality 

HioiM) 



7y {drip) 



< 



H+ 



l{l + l) 



Therefore, in view of (12.71 it suffices to have 

Af2 



-Hio {M) < 



U (Af) pill 



I > 1. 



2l{l + l) 

1 we need to use up all of the good term over T-L'^ which appears in identity 
(12.5). This is not by accident. Indeed, in [3] we will show that even for ^ = 1 we have another 



Note that for I 



conservation law on This law is not implied only by the degeneracy of the horizon (as is the 
case for the conservation law for Z = of Section 12.1 ) but also uses an additional property of the 
metric tensor on This law then implies that one needs to use up all of precisely this good term 
over T-C^. This is something we can do, since we have not used this term in order to estimate other 
integrals. That was possible via successive use of Hardy inequalities. 

12.5. Estimates for the Second Order Derivatives. We can now prove the statement (1) 
of Theorem [3] of Section 4] First note that L ~ ns^ in region A. Therefore, the theorem follows 
from the estimates (12.3) and (12.4), the estimates that we derived in Section [l2.4| (by taking e 



sufficiently small) and the energy identity (12.5). 

Note that the right hand side of the estimate of statement (1) of Theorem [s] is bounded by 
C /j,^ J^H']''i'j:g + J-^^^ Jl^[Ntp]n'^^. This means that equivalently we can apply N as multiplier 
but also as commutator for frequencies I > 1. As we shall see in the companion paper |3^, such 
commutation does not yield the above estimate for I = 0. One can also obtain similar spacetime 
estimates for spatially compact regions which include A by commuting with T and and applying 
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X as a multiplier. Note that we need to commute with in view of the photon sphere. Note also 
that no commutation with the generators of the Lie algebra so(3) is required. 

13. Integrated Weighted Energy Decay 



We are now in position to eliminate the degeneracy of (11.2). 



Proposition 13.0.1. There exists a uniform constant C > which depends on M and Eq such 
that for all solutions ip of the wave equation which are supported on the frequencies I > 1 we have 



Proof. Immediate from the third Hardy inequality, the statement (1) of Theorem |3] and Theorem 
SH □ 



Note that the above proposition shows that in order to obtain this non-degenerate estimate in 
a neighbourhood of Ti.^ one needs to commute the wave equation with dr and thus require higher 
regularity for ■0- This allows us to conclude that trapping takes place along T-L^ . Remark that all 
frequencies exhibit trapped behaviour in contrast to the photon sphere where only high frequencies 
are trapped. This phenomenon is absent in the non-extreme case. 

The statement (2) of Theorem [s] of Section |4] follows from the above proposition and the results 
of Section ini 

14. Uniform Pointwise Boundedness 
It remains to show that all solutions of the wave equation are uniformly bounded. 

Proof of Theorem^ of Section^ We decompose ij} — ipQ A- ip>i and prove that each projection is 
uniformly bounded by a norm that depends only on initial data. Indeed, by applying the following 
Sobolev inequality on the hypersurfaces we have 



(14.1) U>i\\l^^^^) < C[U>i\\h\^^) + ll^>illff^E.) + IV^^ 



where C depends only on Sq- Note also that we use the Sobolev inequality that does not involve 
the L^-norms of zeroth order terms. We observe that the vector field dv — dr is timelike since 
g {dy — dr, — dr) ~ —D — 2 and, therefore, by an elliptic estimate there exists a uniform positive 
constant C which depends on M and Sg such that 



It remains to derive a pointwise bound for ^q. However, from the 1-dimensional Sobolev inequality 
we have 

A^r^{r^,Lo)<- [ J!![M<^. 

ro Js^n{r>ro} 

where C is a constant that depends only on M and Eq. This completes the proof of the uniform 
boundedness of V'- D 

15. Acknowledgements 

I would like to thank Mihalis Dafermos for introducing to me the problem and for his teaching 
and advice. I also thank Igor Rodnianski for sharing useful insights. I am supported by a Bodossaki 
Grant and a grant from the European Research Council. 



LINEAR STABILITY AND INSTABILITY OF EXTREME REISSNER-NORDSTROM I 



35 



Appendix A. Useful Reissner-Nordstrom Computations 



A.l. The Wave Operator. The wave operator in {v 



, I , w , > 



coordinates is 



Dg^p — Ddrdrip + 2dydrip + -dv^ + RdrTp + ^-01 



where i? = £»' + ^ and £>' = ^. 



In (m, v) coordinates 



D' 



Dg^j = - — dud^ir^) V + 

Dr r 

A. 2. The Non-Negativity of the Energy-Momentum Tensor T. We use the coordinate 
system {v,r) and suppose that V = {V" ,0,0) , n = {n" ,0,0). The reader can easily verify 
that if — 2{v^f ^^'^ similarly for n then 

(A.l) 




,9{V.n) 



Appendix B. Stokes' Theorem on Lorentzian Manifolds 

If 72. is a pseudo-Riemannian manifold and P is a vector field on it then we have the divergence 
identity J.^ V^P^ = P ■ nQTz. Both integrals are taken with respect to the induced volume 
form. Note that ng-fi is the unit normal to dTZ and its direction depends on the convention of the 
signature of the metric. For Lorentzian metrics with signature (— , +, +, +) the vector ng-jz is the 
inward directed unit normal to dTZ in case dTZ is spacelike and the outward directed unit normal 
in case 97?, is timelike. If dTZ (or a piece of it) is null then we take a past (future) directed null 
normal to dTZ if it is future (past) boundary. The following diagram is embedded in E^"*"^ 



T 
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